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INTRODUCTION 








There have been a number of studies over the past few years 
regarding the problems associated with coupling two helicopters 
together in such a way as to lift a single load, thus increasing 
the load that can be carried with a given helicopter (References 
1-5). A specific configuration geometry shown in Figure 1 and 
referred to as the twin-lift system is of interest here. 
Preliminary flight demonstrations of this concept described in 
Reference 4 have indicated that while the configuration 
described here can be flown, the pilot workload is high. A 
simplified set of equations of motion were formulated for this 
dynamic system, describing what was assumed to represent the most 
critical aspect of the dynamics and control of such a system. An 
investigation of this reduced system is presented in Reference 1. 
Among other results, it was found that the stability of the 
system is critically dependent upon the location of the tether 
attachment point relative to the center of gravity of the 
helicopters and that a fairly rapidly divergent mode is 
associated with the attachment point located below the center of 
gravity of the helicopter. Further details of this study can be 
found in Reference 1. 

This report extends these studies and presents the complete 
linearized equation of motion for the twin-lift system about a 
hovering trim condition and discusses the features of the 
additional modes of motion involved in the complete system that 
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are not included in the formulation of Reference 1. 


Two sets of equations of motion are presented in this 
report. Part I presents the equations of motion of the twin lift 
system assuming that the yaw angles of the helicopter are 
constant and that the helicopters have a specific orientation 
relative to the load system. Part II generalizes these results 
to permit each helicopter to have a yaw degree of freedom and 
more general geometric characteristics than in Part I. The 
derivation of the equations of motion is described in each 
section and the equations of motion are presented. It is assumed 
that both the helicopters, referred to as the master and slave, 
are identical and that the hovering trim condition consists of a 
rectangular frame orientation, with the tethers vertical as shown 
in Figure 1. Certain simplifications are made in the aerodynamic 
derivatives that are used to describe the helicopter 
aerodynamics. However it should be evident from the formulation 
how to add a more complete aerodynamic description of the 
helicopter to these equations if desired. One of the objectives 
of this study as well as the earlier one described in Reference 1 
was to obtain physical insight into the important features of 
this dynamic system and in particular to note how the system 
differs from dynamics and control problems of a single helicopter 
carrying a single load, a problem that has been exhaustively 
studied in the past. 

Using an aerodynamic description of the helicopters that 
neglects the aerodynamic derivatives that couple the 
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longitudinal and lateral-directional dynamics it can be shown 
that in the case in which the tethers of both helicopters are 
equal in length, the symmetry of the system can be used to 
decouple the equations of motion of this complex system into 
simpler systems which can be used to gain physical insight and 
draw analogues with other systems. In particular it is shown 
that for the linearized system with the assumptions noted above 
and the helicopters oriented initially either with their 
longitudinal axes parallel or perpendicular to the spreader bar 
that the complete equations of motion presented in Part II 
reduced to the following sets: 

The complete sixteen degree of freedom (DOF) system 
described in Part II reduces to four simpler systems: 

1. Yaw motion of Slave Helicopter (1 DOF) 

2. Yaw motion of Master Helicopter (1 DOF) 

3. Planar motion equations (7 DOF). This corresponds 
to system motion in the plane of the paper in 
Figure 1. (Reference 1) 

4. Non-Planar Motion equations (7 DOF). This is 
system motion out of the plane of the paper in 
Figure 1. 

1. and 2. will not be discussed further as they are simple first 
order systems. The planar motion equations are the system 
examined at length in Reference 1. Features of the non-planar 
motion are considered in some detail later in the report. It 
also can be shown that symmetry can be used to divide the planar 


and non-planar motion dynamics into simpler systems. Both the 
planar and non-planar systems can be divided into anti-symmetric 
and symmetric motions. The anti-symmetric motions involve "in- 
phase" motion of both helicopters and are essentially the modes 
involved in maneuvering the system from place-to-place . In the 
planar case this is a f our-degree-of-f reedom system which bears 
some similarities to a single helicopter with a sling load as can 
be seen from Reference 1. Two symmetric systems are present, one 
involving one degree of freedom vertical translation similar to 
the vertical translation of a single helicopter. The other 
symmetric motion involves two degrees of freedom and essentially 
is the same as the dynamics of a helicopter tethered to a fixed 
point on the ground (6). The stability characteristics of this 
motion are very sensitive to the location of the attachment point 
relative to the center of gravity of the helicopter. 

The non-planar motion can also be divided into anti- 
symmetric and symmetric motions. The anti-symmetric motion (four 
degree of freedom) again involves "in-phase" motion of the 
helicopters and translation of the system as a whole, and is 
essentially the same as that of a single helicopter carrying a 
sling load (7). The sling load is a compound pendulum with the 
spreader bar mass and the load mass being the two loads. The 
symmetric motion (three degrees of freedom) involves rotation of 
the entire system about a vertical axis through the load. The 
non-planar motion is considered in more detail later in this 
report . 
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It is interesting also to note that this division into anti- 
symmetric and symmetric modes also indicates a separation of the 
control actions. If the pilots of the master and slave 
helicopter apply exactly equal inputs system motion occurs only 
in the anti-symmetric modes, the system motion is that required 
to move the system from place to place. If exactly equal and 
opposite inputs are applied to the master and slave helicopters 
then only symmetric modes are excited and the total system moves 
about a fixed point in space. 

In general it may be desirable to consider a system in which 
the tether lengths are unequal in which case the decoupling into 
symmetric and anti-symmetric modes does not apply. However 
studies of the unequal tether case in Reference 1 do not indicate 
a strong influence of unequal tether length on the modes of 
motion . 
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NOMENCLATURE 


Coordinates used in the development of the equations of 
motion are shown in Figures 1-4. Generally the following 
subscripts are used to distinguish slave characteristics from 
master characteristics. 

slave associated quantities 


< V < >1 


( ) M . ( ) 2 


( ) T 

< } B 

Geometry of Load 

K ), A( ) 

L 

H 1 

Ho 


h i 
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h’ , H 
Masses 

m h 


master associated quantities 
associated with tether 
associated with spreader bar 

sum and difference coordinates 
spreader bar length, ft 

slave tether length, ft 
master tether length, ft 

displacement of load below spreader bar, ft 

slave attachment point - center of gravity 
spacing, ft 

master attachment point - center of gravity 
space, ft 

are used in the case of equal slave and 
master dimensions, ft 

helicopter mass, slugs 
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load mass, slugs 


M, 


M 


B 


CMB 

Dimensionless 


spreader bar mass, slugs 
spreader bar moment of inertia, 
Parameters 


slug-f t^ 


M = 


m b + m l 


2M 


H 


M. 


P = 


M L + m b 


41 


6 = 


CMB 


(m b + m l )l' 


h i M H 


yi 


Superscripts 

( >* 


( ) 


aerodynamic force derivatives in frame axes 
Moment derivatives with respect to attachment 
point. Velocities of helicopter center of 
gravi ty . 

Stability derivative not divided by mass 

or inertia (X = —) 
u du 


The standard notation for stability derivatives is used, e.g., 



= 1^ ax 

- m h au 

In this notation moment derivatives are relative to the 
helicopter center of gravity. 
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EQUATIONS OF MOTION: PART I 

(Fourteen Degrees of Freedom) 


Since the equations of motion for this dynamic system are 
quite lengthy, the derivation of a reduced set of equations of 
motion is described in this section and then extended to the 
general case in Part II. Here the tether attachment points are 
assumed to be located at the helicopter’s center of gravity and 
the yaw angles of both the master and slave helicopters are 
assumed to be zero. That is, the yaw orientation of both 
helicopters remains constant at an initial orientation in space 
which is as shown in Figure 1. In Part II, both of these 
restrictions on the equations of motion are removed. The 
helicopters may be rotated 90 degrees from the orientation in 
Figure 1 by interchanging aerodynamic stability derivatives as 
noted later (1). Arbitrary initial yaw angles are included in 
the latter formulation. Previous studies of a seven degree of 
freedom model of this dynamic system are described in Reference 
1. This set M degrees of freedom is referred to as the planar 
case. Motion is permitted only in the plane of the paper in 
Figure 1. This section therefore, describes the development of 
the equations of motion including motions out of the plane of 
Figure 1. 

The coordinates used to develop the equations of motion are 
shown in Figure 2. The choice of a relative coordinate system is 
particularly convenient if a Lagrangian formulation is used. A 
space reference frame ( Xg , Yg , Zg ) is defined, with Xg parallel 
to the local gravity vector. The following coordinates are 
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of the slave 


selected 

(x s 


as degrees 
ys , z s ) 


Jn i 0T 1 


<X>B 1 , «1»B 1 


0L 


$T 2 , 0T 2 


of freedom: 

translational displacements 
helicopter attachment point (center of mass) 
with respect to space, zs is parallel to the 
local gravity vector and xs is a displacement 
parallel to the longitudinal axis of the 
slave helicopter. The orientation of the 
longitudinal axis of the slave helicopter 
remains fixed in space. 

Euler angles describing the orientation of 
the slave helicopter tether from the G axis 
system. $ti is a rotation about the Xg axis 
and 0 t i is the second rotation about the 
deflected axis (Yi ' ) 

Euler angles describing the orientation of 
the spreader bar with respect to space. $b i 
is a rotation about an axis parallel to Xg , 
and i is a rotation about the deflected 
axis Zb ' . 

Euler angle describing the orientation of the 
load triangle. 0 l is a rotation about the 
spreader bar axis Yb where Yb is parallel to 
the spreader bar axis. 

Euler angles describing the orientation of 
the tether to the master helicopter. 4 t 2 is 
a rotation about an axis parallel to Xg and 

0 T 2 is a rotation about the deflected axis ip'> 

2 
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These ten degrees of freedom taken with the pitch and roll angles 
of each helicopter 0m, <J)m , 0s, comprise the fourteen degrees 

of freedom of the system. The subscripts M and s refer to the 
master and slave respectively. 

The displacements of the four masses involved in the system, 
the mass of each helicopter Mh , the mass of the spreader bar Mb 
and the mass of the load Ml can be written in terms of the 
geometry of the sling, 

Hi the tether length of the slave helicopter 

H 2 the tether length of the master helicopter 

L the spreader bar length 

Z the distance of the load mass below the spreader bar 
and the rotation matrices 


1 0 0 


0 

cs 

S(|) 


0 

-s <t> 

C(fr 


c 0 

0 

-s 0 

c 

0 

1 

0 

s = 

s 9 

0 

C0 

- 

ctf; 

S(|l 

0 " 


-sij> 

Cl|l 

0 



0 0 1 


The displacements of each mass are 
a.) slave helicopter center of mass 


= cos 
sin 




b.) bottom of slave tether 


< p l> ' < P s> + "’ll *«3 


c.) spreader bar center of mass 


0 


{P CMB } = {P 1> + ^B1 *81^ j^ 2 


d.) bottom of master tether 

(P 2 ) = ( p 2 ) + C*bi ®B1^ 



e.) master helicopter center of mass (attachment point) 

m ( 0 


( P M^ ~ ^ P 2^ + ^ 0 T2 $ T2^ 


0 

H, 


f . ) load center of mass 


{P CML } {P CMB } + [0 L ^B1 ®B1 J 

The kinetic energy of the system due to translation then can 
be written as 

T - \ M h [(P s> * < P s> T ) + \ M H « P M> x {P M )Tl 


+ 2 m b ^ p cmb^ x ^ p cmb^ ^ + 2 m l ^ p cml^ x ^ p cml^ ^ 


CMB 


2 L L 1 CML J 


To this must be added the kinetic energy due to the moments of 
inertia of the helicopters and the moment of inertia of the 
spreader bar. The contributions of the helicopters assuming 
small angles are 
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0 


1 * 2 * 2 1 * 2 '2 
t du = 4- I (e„ + o + 4- i (<t>„ + <t> ) 

RH 2 y M s 2 x M s 


It is assumed that both the master and slave helicopter are 
identical. The contribution of the moment of inertia of the 
spreader bar assuming it to be a slender rod such that Iyy = 
and Ixx = I z 2 = Icmb is 


T RB 2 I CMB ^®B1 COSl ^B 1 ^ + ^Bl ] 


The potential energy of the system is 

V = -Mjjg z - Mgg (z-component of P CM g) 

- M^g (z-component of P^,^) “ MjjS (z-component of P^) 

Lagrange’s equation may now be used to develop the equations of 
motion. Since a large amount of algebra is involved, MACYSMA was 
used to determine the equations of motion. The MACYSMA 
formulation treated only the translational kinetic energy terms. 
The equations of motion were linearized about an arbitrary 
initial geometric configuration of the system assumed that all 
initial velocities are zero, that is only the hovering trim state 
is considered. Certain simplifications in the kinetic energy 
expression arise due to the assumption that all initial 
velocities are zero. In this case, the derivative of the kinetic 
energy with respect to displacements involves only products of 
velocities and thus will be zero in the linearized hover case, 
i . e . , 
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The following names were used for the variables in the 
MACYSMA program 

tt kinetic energy 

faitl, <t>T 1 9 initial value faitlO 

f a i 1 2 , <1>t 2 

f aib i , <J>b i 

saib2, iJjb i 

satatl, 8 t i 

sitat2, 0 T 2 

sitaJL, 0 l 

Other symbols relate directly to the notation above. Thus 
MACYSMA is used to form the kinetic energy and then calculate 
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The potential energy derivatives are also calculated. Helicopter 
body rotation is uncoupled from this procedure and these degrees 
of freedom will give rise to four additional equations. The 
spreader bar moment of inertia gives contributions to the iI»b i and 
$b i equations 


A _ T V 

A dt . T CMB ^B1 


d<t> 


B 1 


A ft } = I CMB *B1 


ao> 


B 1 
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The effect of the external aerodynamic forces acting on the 


helicopter is obtained by calculating the virtual work, 


d , dL ' _ 8L _ .. 

dt • ' aq, ~ k 


where L = T - V 
and 


6 W 


* Q k *qk 


The external forces acting on the helicopter are expressed as 

X , Y , Z , X M , Y m , Z m 
s s s M M M 

Note that these forces must be expressed in gravity axes and not 
body axes. The virtual work of the aerodynamic forces is: 


6W = X 6x + Y <Sy + Z 6z + X M «5x M + Y M 6y M + Z M 6z M 
ss s s ss MM M M MM 


The displacements of the master helicopter must be expressed in 
terms of the coordinates. For small angular displacements: 



II 

<o 

6x 

s 

+ H i 

66 ti 

+ L 

6 %1 - 

H 2 

69 

T2 

• 


6y s 

- H i 

T1 

+ H 2 

6<t) T2 





6z m = 

6z 

s 

- L 

6 *bi - 

' H 1 

[e T1 «e 

T 1 

+ 

*T1 6<|) T1 ] 




+ H 2 

[0 t2 

60 T2 

V T2 

66 

T2 J 



Second order terms are retained in 6zm because the vertical 
aerodynamic forces Zm and Zs have non zero values in the 
equilibrium hovering state. The X and Y forces are zero in 
equilibrium. The equilibrium values of the Z forces are denoted 
Zm o and Zs o . Thus 
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6V = (AX S 4 4X m ) 6 x s 4 (AY m 4 M.) 6y s 4 (2„ 4 Z. ) «z s 

+ ( ' h i ay m - "i Z M 0 *n> «*n + < h i ax m - h i z mo°ii> <e n 

+ (-L Z M ) 64> b1 + (liX M ) 6* B1 4 (0) 69 l 

4 (H 2 AY m 4 H 2 Z mq *„) «» i2 4 (“H 2 &X m 4 H 2 Z mq 8 t2 ) « t2 

The equilibrium values of the vertical aerodynamic forces for the 
rectangular equilibrium state where all initial values of the 
angular coordinates are zero can be seen by inspection to be 

M + M 

Z M0 ‘ - 8 1M H + -V 1 

M L + m b 

Z s0 = - * ' M H + -S -) 

Thus the virtual work contributions to the equations of motion 
are 



X 

s 

iX M + 

AX 

s 

• 

y s 

iV M + 

AY 

s 


z 

s 

Z M0 * 

Z s0 + 

• 

®T1 

- h i (ay m 

+ Z M0 


0 T1 

Hj(M m 

Z M0 

• 

$ B1 

- l(Z M0 

+ az m ) 


*B1 

liX M 


• 

0 L 

0 





H 


'M 


M L + m b. 


m l + m b 


T1 


’T1 


m l + m b 


*H 


) - LAZ 


M 
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- gH 2 (M H + 


<j> T2 H 2^ AY M + Z M0 


M L + M B 


> *T2 + H 2 AY M 


0 T2 H 2^ AX M Z M0 0 T2^ ~ gH 2^ M H + 2 ^ 0 T2 H 2 AX M 

Note that some constant terms appear which will be balanced by 
terms on the right hand side and in addition there are 
contributions to the spring matrix due to the initial value of 
the Z-force. Only the dominant aerodynamic derivatives are 
included in this mode. The perturbations in the aerodynamic 
forces in a gravity fixed reference frame are: 


AX = X u + X nl B, - T. 9 
u B1 1 0 


A y - Y y v + Y a1 Al + T 0 * 


AZ = Z w + Z Q 9 
w 0C c 


For the slave helicopter 


u = x 
s s 


= y. 


W = 2 

S S 


For the master helicopter 


U M X s H 2 0 T2 + H 1 0 T1 + L ***B 1 
V M = Y s + H 2 ®T2 “ H 1 ®T2 


w w = z - L G> , 
M s B1 
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The four equations describing the helicopter body motions are 

I 9y = AM m 
y M M 


I <t> = AL m 

x M M 


I 8 = AM 

y s s 


where 


I 

X 

M s 




AM 

= M u + 
u 

M q + 

q 

S B1 

B 1 

AL 

= L v + 

V 

L p + 
P 

^1 

A 1 


where for simplicity aerodynamic coupling derivatives have been 
neglected . 

Thus using these results with the output of the MACYSMA 
program the equations of motion can be written as in general as, 

[M] {x} + [C] {*} + K {x} = B {u} 

where the fourteen motion variables are: 

{x} = U s .y s «z s .* T1 . 0 T1 » < t» B1 «^Bl’ 0 L ,<t> T2 ,0 T2’ 0 M ,0 s ,<,, M ,<,> s ] 
and the controls are: 

{u} = (B ls > b 1m , A lg , a 1m , e Cg , e CM ] 

The [M] , [C], [K] and [B] matrices can be obtained for this 

case from the more general matrices given in the next section. 

The equations have been normalized by mass and inertia. With the 
exception of the 0 l equation, the first ten equations have been 
divided by the helicopter mass (Mh ) . The load coordinate ( 0 l ) 


equation is divided by the load mass (Ml), The pitch and roll 
moment equations are divided by the respective moments of inerti* 
of the helicopter. 

The following nondimensional quantities are introduced: 


m b + m l 

p = — > ratio of spreader bar plus load mass 


H 


to total helicopter mass 


Note that 


P = 


M t 


ratio of spreader bar mass to load 


+ mass plus spreader bar mass. 


M, 


(1 - P) = 


M L + m b 


M, 


(P) (1 - P) = 


2M 


H 


It can be noted from the M, C, and K matrices that the equations 
of motion for the complete fourteen degr ee-of-f reedom system 
decouple into two seven degree-of-f reedom systems, the planar 
dynamics studied in Reference 1 associated with the following 
degrees of freedom: 

* X p* = ^ y s’ Z s’ ^Tl’ ^B 1 ’ ^T2 ’ *M’ *8^ 

The remaining coordinates describe what will be referred to as 
the non-planar motion. It is associated with the following 


degr ees of freedom 


Note that a slightly different notation is used in Reference 1 
for the coordinates. 



EQUATIONS OF MOTION: PART II 


(Sixteen Degrees of Freedom) 

The previous section described in detail the derivation of 
the equations of motion for the twin lift system. The 
formulation in Part I is restricted to the tether attachment 
points of both helicopters located at their respective centers of 
gravity. In addition the initial yaw orientation of each 
helicopter is assumed to be either perpendicular or parallel to 
the load suspension frame and the yaw angles are assumed to 
remain zero throughout the ensuing perturbed motion. The 
formulation is given for in Part I for the longitudinal axes of 
the helicopters perpendicular to the spreader bar in the initial 
state. It can be seen from Part I that the case in which the 
helicopters are rotated 90* is readily treated by interchanging 
the aerodynamic stability derivatives. 

In this section, the formulation is extended to include 
vertical spacing between the tether attachment points and the 
centers of gravity of each helicopter. This dimension was shown 
to have an important effect on the stability of the twin lift 
system in the planar case. In addition, the initial orientation 
of the master and slave helicopters with respect to the load 
frame is arbitrary as shown in Figure 3. The following 
quantities are added in this section: 

4>s o trim or initial yaw angle of slave helicopter relative 
to frame (spreader bar) 

i|>Mo trim or initial yaw angle of master helicopter relative 
to frame (spreader bar) 




A4)s change in yaw angle of slave helicopter measured from 
initial orientation of slave helicopter. 

Ai|jm change in yaw angle of master helicopter relative to 
initial orientation of master helicopter. 

The equations of motion are formulated such that it is assumed 
that the yaw motion of the master relative to the slave 
helicopter , 


Aifi 


MR 


A *M 


M's 


is a small quantity. The change in yaw angle of the slave 
helicopter is not restricted to be small. 

Two equations of motion are added to the dynamic system 
described in Part I. 


A* s - N r A*. = N « IR 

in S 


and 


A*M - N r A* m = N 6 

IK M 


These modifications to the equations of motion can be 
incorporated into the equations of motion in the following way. 
A frame axis system (Xf, Yf , Zf ) is added which moves about Zf 
with the change in yaw angle of the slave helicopter as shown f 

Figure 4. Within the framework of the small angle assumption' 

/ 

made in the linearization of the equations of motion, 4 >b i ca/ 
interpreted as the rotational displacement of the spreader 
relative to space. Note however, that the other angular 
coordinates describing the system motion rotate with Adis 
helicopter body axes are now misaligned from the frame e 



Figure 4, the relative spreader bar motion is also shown. Thus 
to summarize: 

(Xf, Yf , Zf ) Frame axes rotate with Atjjs . The frame motion 
variables defined in Part I are perturbations 
relative to this axis system. (Figure 4) 

(Xbs, Ybs, Zbs ) Slave helicopter body axes. Misaligned from 

frame axes by tfis o and rotate with AtJjs and 
thus remain in a fixed orientation relative 
to the frame axes. 

(Xbm, Ybm, Zbm) Master helicopter body axes. Misaligned from 

frame axes by i|iho and rotate with A<J;m . Thus 
it is assumed that in the disturbed motion 
that (AiJjm - A»Jjs ) is a small quantity. 

The equations of motion are now modified to account for the 
misalignment of the helicopter body axes and the frame axes. The 
stability derivatives of the helicopter are given in terms of 
velocities and forces oriented in the body axis directions. Note 
that these are not true body axes as they do not rotate with 
either the pitch or roll angles and thus the thrust vector must 
be accounted for in the expressions for the longitudinal and 
lateral aerodynamic forces. In addition, the equations must 
account for the fact that the tether attachment point is being 
used as the reference frame for the translational velocities. 

First, the transformation of the aerodynamic derivatives is 
considered and then additional inertial and gravity terms due to 
the spacing between the attachment point and the center of 
gravity of each helicopter are considered. In this formulation 
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only the dominant aerodynamic derivatives of the helicopter are 
considered. The formulation indicates how additional terms may 
be included. It is assumed that the following aerodynamic 
derivatives characterize the helicopter: 



The following control derivatives are included 

X B1 ' M B1 * Y A1 ’ L A1 ’ N 6 * Z 0c 

I K 

The velocities of the helicopter center of gravity expressed in 
terms of attachment point velocities (ut, vt ) and attachment 
point, center of gravity spacing (h’) are: 


u = u T - h’ 6 
v = v T + h’ <t> 

The final formulation permits the spacing between the attachment 
point and the center of gravity to be different on each 
helicopter with the subscript 1 referring to the slave and the 
subscript 2 referring to the master. 

The aerodynamic forces in the helicopter body/gravity axes 

are : 


AX_ = X u + X nl B, - T 8 
B u B1 1 o 


AY = Y v + Y.. A, + T <t> 
B v A1 1 o 


The body motion equations are about the helicopter center of 
gravity are: 


I 8=Mu+M8+M n1 B, i 
y u q B 1 1 
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V = V + V * l n A 1 

The moment equations will be retained in body axes and the 
force equation will be transferred to frame axes. 


ax f = 

ax b 

cosi|; 

- AYg sini|> 

> 

*< 

II 

ay b 

COSlJj 

+ AX D s i n iJj 

D 

U T 

X 

s 

cosiJj + 

y sinijj 
s 

V T 

X s 

COSlli “ 

x s in»Jj 
s 


In the force equations, the effect of the spacing between 
the attachment point and the center of gravity can be neglected. 
This is consistant with neglecting the stability derivatives Xq 
and Y P . 

Thus using the above transformation 

AX F = X u *S + X V y s + ( X Bi cos ^) B ! - (Y A1 sinOO A l 

+ T o (-8 cos ijj - sinJO 

AYp = Y* x g + Y v *y g + (XgjSiniJj) Bj + (Y^cosiji) Aj 

+ T q (~8 sinil» + <b cos ip) 

where 

X* = X cos^i|j + Y sin^ijj 
u u v 

X* = (X u - Y v ) sinijj costjj 

Y* = (X - Y ) simj> cosi)> 
u u v 

Y* = Y v cos^'Jj + X y sin^ijj 
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If it is assumed that X = Y then 

u v 



X* = Y* = 0 
v u 

This is a reasonable simplification and is used in equations 
of motion presented later. 

The linearized form of the aerodynamic forces is obtained 
by letting ^ = i|Js o for the slave helicopter and i|j = i|>mo for the 
master helicopter terms. 

To account for the spacing between the attachment point and 
the helicopter center of gravity, the moment equations are 
written about the attachment point. Thus the stability 
derivatives about the attachment point are denoted 

M* = M - h’ AXg ’ 

L* = L + h’ AYg ’ 

Note that AX * and AY * are body axis forces and do not contain 
B B 

the thrust vector. 

The moment equations for the helicopter written about the 
attachment point are: 


[I + 
y 

h’ 2 

V 

«• 

0 - 

h MgU 

- Wh’ 

e = M 

[i + 

X 

h’ 2 

M h ) 

»• 

4> + 

h ’v 

- Wh’ 

d> = L 


The following terms are added to the equations of motion to 
account for the displacement of the helicopter center of gravity 



from the attachment point. The inertial terms are: 


AX 

(-=*> 

m h 


-h„ [cosik, 0 M + sinik, ] -h. [cos»J> 0 + siniji <|> ] 

2 Mo M Mo M 1 so s so s 


AY 

(— ) 
V M J 

H 


-h 0 [cosik. 4).. - sirnk.. 0..] +h. [costf/ <t> - siml) 0 ] 

2 Mo M Mo M 1 so s so s 


^T 1 

(— LL) 

'-M 7 

H 


h 2 H l ( - cos ^Mo <t, M + sinil, Mo e M ) 


8 T1 

( M 1I) 

H 


h 2 H i ^ -cos ^Mo^M " si "*MoV 


(il) 

(m H 


" h 2 L ( cos ^Mo 8 M + sinq, Mo*M ) 


*^T 2 

V M ; 


h 2 H 2 < co ‘*Mo*M - sin *Mo*M ) 


8 T2 

(r r^) 

m h 


h 2 H 2 (cos,1, Mo e M + s ^ nl ^Mo V 


y 


(h 2 r> <-*« COSl * , Mo * !i "*«o + H l ( ‘®n cos,|, Mo 

y 

+ *T1 sina; Mo ) + V®T2 cos *Mo ~ *^T2 sin %o ) 


+ L(-4» fil cos 4» m<j ) + h’0 M } 


(?) 

X 


M 

<h 2 ?> + *s cos,l, Mo + "l'Al sin *Mo 

X 


~*T1 cosa, Mo } + H 2 (8 T2 s in ^Mo + ^T2 cos, W 


+ L(-4» b1 sin* Mo ) + h* <D m } 


0 


h l ? { -*s cos *so - si,, Ko 4 h l ®s } 

y 


4> 


m h •• 

h, t — {-x sinjj + y cosiJj + h, <J) } 
II 1 s so J s so Is 

x 


In addition the gravity terms must be resolved as shown in 
the equations of motion and the following terms added to the body 
moment equations: 


0 


M 


M H 

g h 2 T~ 0 M 


4>, 


M 


M 


H 


g h 2 I '•'M 
x 


<k 


e 


m h 

- 8 h i - e s 
y 


i 


m h 

-« h ir*s 

x 


Introducing these terms and the resolution effects for helicopter 
orientation as noted leads to equations of motion or the form 

[M] x + [Cj x + [K] x = [B] u 

where 

x = <x s , y g , z s , 4* t1 , 0 t1> 0 L . 4 > t2 i 0 T 2 * e M’ 0 s 

V V V V* 

T 

U = {B ls’ B 1M’ A 1 s* A 1M’ 0 cs’ 9 cM’ 6 TRM’ 6 TRs* 


The system is described by sixteen degrees of freedom and there 
are eight controls. The matrix elements are given on the 
following pages. The yaw moment equations are gien in the text 
and not included in the matrices. The equations of motion for 
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the simpler case described in Part I are obtained by choosing Jjs o 
= i|)mo = 0, or i|jso - = tt/2, and eliminating the yaw moment 


equat ions . 
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NON-PLANAR MOTION 


It can be seen from the equations of motion developed in the 
previous sections that the complete equations of motion for the 
twin lift system decouple into independent sets of equations of 
motion if the lateral-longitudinal coupling terms due to 
aerodynamic stability derivatives are neglected in the 
helicopter models and the initial yaw angle of each helicopter is 
either 0 or 90" (longitudinal axis of the helicopters 
perpendicular or parallel to the spreader bar). 

The uncoupled sets of equations of motion involve the 
following degrees-of-f reedom if the initial orientation of the 
helicopters is such that their longitudinal axes are 
perpendicular to the frame (iJjmo = ifiso = 0). 

Planar Motion (seven degrees of freedom) 

{x P } = {y s> V *T1’ i * < t > T2 > <V V T 

Non-planar Motion (seven degrees of freedom) 

^ X Np ^ = {x s> 0 T1’ 0 T2’ 1 ’ 0 L’ 0 M’ 0 s* 

The other two degrees of freedom, the yaw angles of each 
helicopter are uncoupled. If the initial yaw angles are 90", / 

then the pitch and roll angles are interchanged in the degrees 
freedom of the two cases. 

The planar motion has been investigated extensively and 
results of this investigation are presented in Reference 1. 
this section the basic dynamics of the non-planar case are 
examined. The same physical parameters are used as in Re 



I 


1. These are listed in Table I. The helicopters used in the 
example correspond nominally to the UH-60, Blackhawk. 

It was noted in the studies of the planar motion that if the 
tether lengths are equal (H = Hi = H 2 ) and the attachment point 
center of gravity spacing is the same on both helicopters, then 
the dynamic system can be further subdivided into sets of 
symmetric and anti-symmetric motions and considerable insight 
gained regarding the dynamics of this complex system. The anti- 
symmetric modes involve the two helicopters moving "in-phase" and 
are essentially the modes involved in translating the entire 
system laterally, i.e., in the plane of the paper in Figure 1. 

The symmetric modes involve essentially "out-of-phase" motion of 
the two helicopters with the load stationary and are equivalent 
to those of a single helicopter tethered to a fixed point. See 
Reference 1 for additional discussion. 

A similar division or decoupling of the modes can be made in 
the non-planar case. Again using the terminology anti-symmetric/ 
modes to refer to the case in wh 
phase", the following coordinate 

Anti-symmetric motion: 

*s’ 0 T1 = 0 T2’ 0 L’ 0 M 

This is a four degree of freedom 
are those of a helicopter with a 
compound pendulum (7). Figure 5 
motion . 


ich the two helicopters move "i; 


relationships are involved: 


j 

/ 


V <*81 * °> 


system. The equations o 
sling load that acts li 
shows the geometry of 
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Symmetric motion: 


Oj 


B 1 


T 1 


= - 9 


T2 


M 


= - 0 


(e T = 0) 


This is a three degree of freedom motion that involves only 
rotation of the system about a vertical axis through the load, 
spreader bar center of mass as si:- -v;n in Figure 6. 

The next section shows how the equations are decoupled by 
the introduction of new coordinates. 

Decoupling the Equations of Motion (Non-Planar case) 

The equations of motion for the non-planar case can be 
written as 

[M] { x N p} + [C] {* Np } + [K] {* Np } = [ B ] {u} (1) 

where 

{* N p} = [x g , ® T1 » e T 2> 1 ’ 0 L’ e M’ 9 S ] 

The matrices [M], [C], [K] and [B] are given following this 
section for this reduced degree of freedom case. 

The case in which the tethers are equal (H = Hi = H 2 ) is 
considered and a new set of coordinates similar to those used in 
the planar case is introduced to decouple the equations. The 
decoupling procedure is as follows. 

Noting that the translational displacement of the master 
helicopter can be written as 


X M X s + H(9 T1 9 T2 ) + LJj B1 


(3) 


This coordinate is introduced and sum and difference coordinates 


are defined: 


Ix = 


ze = 


20 , 


X 

s 

+ 


Ax = 


2 


0 

S 

+ 


A0 = 


2 


0 T1 

+ 

0 T2 

A0 T = 


2 



* • m O 

M s 


M s 


( 4 ) 


The load displacement 0 l is chosen as the seventh coordinate. 
The spreader bar yaw motion tiiB 1 can be eliminated as a motion 
variable, using the equation for xm given in equations (3), 


. Ax A H .. 

*B1 r + L i9 T 


Sum coordinates are associated with the anti-symmetric motion and 
difference coordinates are associated with the symmetric motion. 
The pitching moment equations for the master and slave 
helicopters can be written as 


• t * * * ' I Mil O My (, My X 

®M ~ M u X M “ M q 6 M h T~ X M + h ’ T~ ®M ' gh ’ T~ ®M M B1®1M 

y y y 

(5) 

e - Mx - M*0 - h’ / x + h’^ 0 — gh’ =4 0 c = M*B, 

susqs Is Is Is B1 is 

y y y 


Sum and difference equations are obtained by adding and sub- 
tracting these equations. 

u ± ± * ' My i* ? My •* My * 

10 - M Zx - M 10 -h’ ■=?- Ix + h*‘ y 1 I© -gh’ ^4 10 = M* IB. 

U q I I I Bll 

H y y y 

( 6 ) 
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w * * 

A0 - M Ax 
u 


* M u 9 " 

M A0 -h’ ^ Ax + h’* A8 

q I I 

y y 


m h * 

- gh’ f- A9 = M^ABj 
y 


(7) 


Note that the pitching moment derivatives are expressed with 
respect to the tether point. In order to reform the other five 
equations it is convenient to introduce the following set of 
additional coordinates which describe the linearized 
displacements of the centers of mass of the spreader bar and 
load, and the master helicopter tether point. 


X M 

X 

s 

- H 

A0 t 

+ 

L *B1 ’ 

* 

X M = 

X M - h> 

9 M 


X B 

X 

s 

+ H 

9 ti 

+ 

L . 

2 *^B 1 ’ 

* 

X = 

s 

x - h’ 
s 

9 

s 

(8) 

X L = 

X 

s 

+ H 

9 ti 

+ 

L , 

2 ^Bl 

Ze L 





The equations of motion in terms of these coordinates are: 

x equation, 
s 


M 'x + M u x„ + M *x n + M t x t - X (ft* + x*) + T 219 - h’ M 2Z9 
Hs HM BB LL u s M o H 

(9) 

= *B1 2ZB 1 

9^^ equation divided by H 

M H X M + M B X B + M L X L “ *u X M + g( 2 } 9 T1 + T o 9 M “ h ’ M H ®M 


( 10 ) 


= X B1 B 1M 


equation divided by L 
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( 11 ) 


• • 1 o •• »• 3k •• 

M H X M + 2 (M B X B + M L X L ) + ~L ^Bl ” X U X M + T o 9 M “ h * M “ 9 » 


‘H ~M 


x bi b im 


0 g equation divided by ZM g 

x g + g 0 g = 0 (12) 

©12 equation divided by H 

M + M 

_M H X M + X u X M + g( 2 ^ 6 T2 T o 0 M + h ’ M H 9 M = _X B1 B 1M 


Noting that 


B1 " L ix + L i9 T 


X B ' X L ' ' Z0 l 

Adding the 0 T1 and 0 T2 equations (10 + 13) 

Vb + M L X L + g(M L + M B ) E9 T = ° 

Adding the ii) g ^ and © T2 equations and multiplying by 2, 


(13) 


(14) 


(15) 
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CMB 


Vb + M L*L + -f^ *B1 + g(M L + V 9 T 2 ' 0 


(16) 


It may be noted upon comparing equations (15) amd (16) that if 
the moment of inertia of the spreader bar about its center of 
mass is neglected then 

0 T2 = ^ 0 T^ 

implying that 0^^ = 9^ ^ > that is, the two tether angles two 
are equal if the spreader bar inertia is neglected. In this 
case, one motion variable becomes redundant as does one equation 
of motion. xb and xl can be expressed as from (4), (8) and (14) 


as : - 
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Xg B Sx + H ie T 


Xg = Zx + H E0^, + Z0g 


( 17 ) 


The sura coordinate equations (anti-symmetric equations) are 
formed from the Xs equation (9), and the load equation (12), the 
sura pitching moment equations (6) and the (0ti + 0T2) equation 
(15) • Note that these substitutions transfer the pitching moment 
derivatives to the center of gravity. The anti-symmetric system 
equations of motion are: 


• 

•* # * 

Ex - X Ex - 
u 

h’Z0 + g(l 

+ p) 10 - 

gp£0 T 

= 

x bi zb i 


- M Ex* + 

u 

10 - M EG 
q 

+ e y gP(20 

- Z0 T ) 


= M B1 ZB 1 


It 


H 

II 


( 18 ) 


Zx 


+ HI0 t 

+ Z0 T + 
L 

g0 L 

= 0 

• 



- gS0 T 

it 

+ PZ0g + 

g0 L 

= 0 




Zx 

- h*ze - 

Zx* 

= 0 


The difference 

coordinate 

equat ions 

( symmetric 

mode 

equations) 


are obtained by subtracting two times the \\)b i equation from the 
Xs equation, the difference moment equation, and a combination of 
the 9 t i , 8t2, and iJjb i equations. The symmetric system equations 

of motion are: 


Ax - X^Ax* + h’A8 + g(l + p) A8 - gp AQ^ 


- M Ax* + A0 - M A0 + e gp(A0 - A0_) 
u q y T 


= x bi ab i 


= M B1 AB 1 


•I 

6 Ax 


+ 6HA0 t + g AC‘ t 


Ax - h ’ A0 - Ax 


= 0 
= 0 


( 19 ) 


Thus the new coordinates invoved in these uncoupled sets of 


39 



I 


equations of motion are: 
X 


NP, AS 


{£x, Z0, 10^, 0^} 


X Nps = {Ax, A0, A0 t } 1 

The parameter 6 is a non-dimensional expression of the spreader 
bar moment of inertia. If the spreader bar moment of inertia is 
equal to zero, 6=0, and the last equation in the symmetric 
case gives A0t = 0 and the symmetric set reduces to a two degree 
of freedom system. 

Note that x refers to the displacement of the tether 
attachment point while x* refers to the helicopter center of mass 
displacement. The characteristics of these two systems are now 
considered. 

Anti-symmetric Modes 

It can be noted that the anti-symmetric equations (18) are 
the equations of motion of a helicopter coupled to a sling load 
which is a compound pendulum. The parameter 3 represents the 
ratio of the two supported masses, the spreader bar mass and the 
load mass. In the limit of no spreader bar mass (3 = 0), £0t = 

0l , the load equations of motion reduce to those of a single 
pendulum of length (H + Z). Comparison with the planar anti- 
symmetric equations of motion shows that for h* = 0, Z = 0 (load 
on spreader bar) and 3=0, the equations of motion are identical 
if the following interpretation is made. 

- H(I0 t - Z0) = H(0 l - 10) 

The coupled frequencies of the double pendulum system can be 
determined from the characteristic equation of the two degree of 


40 



The characteristic equation of this 


freedom system for E9 t , 0l . 
system is, 


HZ P 
g(H + 

When 3=0 



and when P 


4 2 

zy s + s + 

(no spreader 

£ 

H + Z 

1 (no load 


g = 

(H + Z) 

bar mass) 


0 

there is a single frequency 


mass) there are two frequencies 



£ 

H 


u 


Figure 7 shows the variat 
p. At the nominal value 
frequencies are: 

= .826 rad/sec 

u g 2 = 8.16 rad/sec 


= £ 

2 Z 

ion of these uncoupled frequencies 
of p, with load (p = .054) the two 


with 


There is a wide variation in one of the frequencies associated 
with the load as the load is reduced to zero (P = 1.0). 

The anti-symmetric equations of motion are those primarily 
associated with maneuvering the system from place to place. The 
complete system involves these two sling modes coupled to the 
helicopter dynamics. For the nominal load (P = .054) the high 
frequency ( cj s 2 ) associated with the sling is only weakly coupled. 
Symmetric Modes 

The symmetric motion as described by equations (19) show as 
noted above that if the spreader bar inertia distribution is 
neglected then 6 -» 0 and A0 t = 0 indicating that the tethers 
remain vertical throughout the motion in this limit. The tether 
cable tension provides a restoring moment as indicated by the 


equations of motion. With the spreader bar inertia included 
there is an angular oscillation associated with the spreader bar. 
The frequency of this motion uncoupled from the helicopters is 
obtained from the last of equations (19): 

2 _ &_ _ U 


u 


<1> H<5 Hp 


where 
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x _ CMB = B 

° 2 3 

(M b + M L )L 2 

The nominal value of P with load is .054 and therefore with H = 

13.25 ft, this frequency is: 

o . = 12.07 rad/sec 
i|J 

The spring in this motion is tension in each tether. Since this 
frequency is high compared to the helicopter dynamics it is 
probably quite a good approximation to determine the influence of 
the spreader bar inertia on the modes associated with the 
helicopter dynamics by using a quasistatic approximation. 


Ax £ " 6 A0 T = " Hu O. 2a0 T 


If the entire system is moving at a frequency w, then the 
differential tether angle is 


( 20 ) 


Aa = ,u ,2 Ax 

A0 t - ( zr } ir 

ijj 


Thus the tethers tend to remain relatively near vertical as 
indicated by the limit 6 ■* 0, and the spreader bar follows the 
helicopter motion, i.e., 



The equations of motion assuming quasi-static bar motion are (P 
<< 1) obtained from equations (19) and (20), 


(1 + 6p ) Ax - X u Ax* - h’A0 + g(l + p) A0 s X B1 AB ^ (21) 

e <5pAx - M Ax* + A0 - M A0 + e gp A0 s M D , AB, 
y u q y B 1 1 

Thus fully loaded with 6 = .018, p = .43 the effect quasi-static 
tether rotation, proportional to 6p, on the helicopter modes is 
quite small. With no load, p is very small and it would also be 
expected that the coupling between the angular motion of the 
spreader bar and the helicopter motion is weak. Note that the 
uncoupled frequency of the spreader bar motion varies 
significantly with load due to the tension variation as shown in 
Figure 8. 

The limiting mode shapes in the symmetric motion with 
variation in spreader bar inertia are shown in Figure 9. It also 
can be seen that in the limit of very large spreader bar inertia, 
Ax = H A0 t and the anti-symmetric equations become identical to 
the symmetric planar case. That is, the spreader bar tends to 
remain fixed in space and a divergent mode will occur for the 
attachment point located below the center of gravity of the 
helicopter. This limit can be readily seen by changing 
coordinates in the difference equations replacing A0 t by i|»b i . 
However the actual system will tend to be far from this limit for 
the physical parameters typical of this system. 
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In the practical case where the non-dimensional parameter 
measuring the spreader bar inertia (6) is small the character- 
istic equation for the anti-symmetric coupled motion, neglecting 
6p terms in equation (21), which is equivalent to assuming that 
the tethers remain vertical is, 

(s) (s - X u ) (s - M q ) + g(l + p) M u + e y gp(s - X y ) = 0 

For e y positive, there is a stabilizing tendency, i.e., since in 
this limit the tethers are vertical an effective attitude 
stability is provided by tether tension acting at the attachment 
point below the center of gravity of the helicopter. Placing the 
attachment point above the center of gravity acts like an 
unstable attitude feedback. This stabilizing trend is opposite 
to that shown for the planar symmetric case in Reference 1. 

Figure 10 shows numerical values for the modes of motion for 
the non-planar system with load (the symmetric and anti-symmetric 
modes) as a function of tether attachment point location relative 
to the helicopter center of mass. The uncoupled helicopter modes 
for this sample calculation based on the parameters in Table I 
are -.82 and +.172 ± .572i. The destabilizing effect of locating 
the attachment point above the center of gravity can be seen from 
the figure. The movement of the complex pair is very similar to 
that obtained with attitude feedback. 
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TABLE I 

NUMERICAL VALUES OF SYSTEM PARAMETERS 


Helicopter Mass Mg = 435 slugs 


Load Mass 
Spreader Bar Mass 
Slave Tether 
Master Tether 


Mjj = 353 slugs 


Mg = 20 slugs 


Hj = 13.25 ft 


H 2 - 13.25 ft 


Spreader Bar Length L = 68.9 ft 

Load distance below spreader bar Z = 34.5 ft 

Tether Attachment point/CG spacing h’ = 3.5 ft (nominal) 

Helicopter moment of inertia in pitch I y = 43,000 slug ft' 

Helicopter Stability Derivatives (referenced to CG) 

x = -.0602 sec * 
u 


M = -.415 sec 

q 


-1 


M = .00538 (ft-sec) 
u 


-1 


Non-dimensional parameters: 


m b + m l 

P = ~4r. = 0.43 


2M 


P = 


6 = 


M L + m b 


054 
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M, 


(M b + M t ) L 


3(M b + M l ) 


= .018 


h’M 


(uniform mass distribution) 


H 


030 
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SUMMARY 


A full set of equations of motion for the twin lift system, 
linearized about a hover trim condition have been derived and 
presented. It is shown that this full set of equations of motion 
decouples into simpler sets of equations of motion if the 
aerodynamic coupling derivatives of the helicopters are 
neglected. One of these decoupled sets of equations of motion 
(referred to as the planar set) was studied at length in 
Reference 1. The other decoupled set (referred to as the non- 
planar set) is examined here. It is shown that when the 
geometric configuration of the twin-lift is symmetric that a 
further decoupling is possible into anti-symmetric and symmetric 
sets of equations. One set of these reduced equations of motion 
referred to as the anti-symmetric set is directly equivalent to 
the longitudinal motion of a single helicopter with a sling load. 
The second set, referred to as the symmetric set corresponds to 
rotation of the entire system without load motion. As in the 
case of the planar symmetric motion, the location of the tether 
attachment point influences the stability of the non-planar 
symmetric mode. The trend is opposite however in the non-planar 
case in that an attachment point below the helicopter center of 
gravity gives a favorable effect on the stability. The effect is 
not as strong however as the unfavorable effect on the planar 
symmetric mode shown in Reference 1. 
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